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Abstract The present work details the analysis of the aero- 
dynamics of an experimental swirl stabilized burner repre- 
sentative of gas turbine combustors. This analysis is carried 
out using High Frequency PIV (HFPIV) measurements in a 
reactive situation. While this information is usually available 
at a rather low rate, temporally resolved PIV measurements 
are necessary to better understand highly turbulent swirled 
flows, which are unsteady by nature. Thanks to recent tech- 
nical improvements, a PIV system working at 12 kHz has 
been developed to study this experimental combustor flow 
field. Statistical quantities of the burner are first obtained 
and analyzed, and the measurement quality is checked, then 
a temporal analysis of the velocity field is carried out, indi- 
cating that large coherent structures periodically appear in 
the combustion chamber. The frequency of these structures 
is very close to the quarter wave mode of the chamber, giv- 
ing a possible explanation for combustion instability cou- 
pling. The frequency of these structures is very close to the 
quarter wave mode of the chamber, giving a possible expla- 
nation for combustion instabilities coupling. 



1 Introduction 

Stronger regulations concerning pollutant emissions led to 
the development of new generations of swirl burners stabi- 
lizing lean premixed flames well suited to reach low NOx 
levels (ICorreaL 119981) . Swirling motion permits to improve 
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the mixing rate between fuel and oxidant streams and to con- 
trol the flame stabilization through the swirl-induced recir- 
culation of hot pro ducts near the nozzle (IChen and DriscollL 
1988HSyredil2006l) . Above a critical swirl number (S > 0.6), 
the expansion of the swirling flow leads to a vortex break- 
down, which is accompanied by an inner recirculation area. 
This zone provides the major mechanism for flame stabiliza- 
tion since this recirculation area of hot combustion products 
continuously supplies energy for the ignition of the incom- 
ing fuel-air mixture. 

However, under lean conditions, these flames tend to ex- 
hibit undesired combustion instabilities, which may cause 
flame extinction, flashback, reduc e engine life time or lead 
to cat astrophic structural damages (ICandelll2002tlNauert et al. 
20071), despite the stabilizing effect of their recirculation zones 
The mechanisms of these instabilities are based on various 
complex interactions between combustor geometry, flow field, 
pressure, mixing, chemical reactions and heat release. It has 
been shown that swirled flows develop characteristic peri- 
odic large coherent structures, which play a n essential role in 
the dy namics of turbulent swirled flames (iPaschereit et al.L 
1998h . Their interactions with the heat release process and 
acoustic resonant modes of the combustor can cause unde- 
sirable thermo-acoustic instabilities in the combustion sys- 
tem. Due to the importance of these structures as drivers of 
combustion instabilities, numerical and experimental stud- 
ies have been undertaken in order to understand their in- 



teractions with chemi stry and acoustics (iRoux et al.L 12005 : 
Sengissen et alll2007l) . 

Lab scale swirl stabilized combustors have been developed 
and studied using microphones, OH* or CH* emission sen- 
sors ([Paschereit et a l., 1999), photographi c imaging techniqu e 



dBroda et al. 



Meier et al 



20061) . 



19981) or laser diagnosti cs ( Janus et all l2005l 



2006l:IOlivani et al.Ll2007l:IWeigand et al.Ll2005l 



Laser-based tools offer the potential to measure most of the 



important quantities necessary to gain a deeper insight into 
complex chemical and physical processes in turbulent flames. 
The flow field can be measured by Laser Doppler Velocime- 
try (LDV) or Particle Imaging Velocimetry (PIV) (IJi and Gore , 
2002t iMueller et all [l998l: IWoodmansee et all boOTJ ). mix- 
ing process and flame structure by Planar Laser-Induced Flu- 
orescence (PLIF), while laser Raman scattering permits to 
determine the major species concentrations, temperature, and 
mixture fraction. Up to now, due to laser technology limita- 
tions, these information were usually sampled with a rep- 
etition rate of an average of 10 Hz (except for LDV), giv- 
ing combustion statistical quantities that significantly helped 
in understanding turbulent flames. However, some of the 
mechanisms leading to combustion instabilities are still not 
well understood, as the interactions existing between large 
coherent structures, chemical reactions and the acoustics of 
the burner are complex and usually strongly dependent on 
the configuration. Indeed the number of experimental stud- 
ies investigating the temporal evolution of the flow field in 
gas turbine combustion chambers is rather limited. This is 
partly due to technological limitations of lasers and cameras, 
which do not permit to acquire flow field data at a high ac- 
quisition rate. However, combustion instabilities, flame ex- 
tinctions or flashbacks are by nature very unsteady phenom- 
ena, whose comprehension is built on the study of spatially 
and temporally resolved data. 

In the present work, the flow field of a lean swirl stabilized 
flame in a gas turbine model combustor is investigated using 
time resolved PIV. For the last 20 years, the development 
and improvement of PIV have proceeded in several stages, 
related to the technical progresses successively achieved in 
the fields of lasers (double-pulsed solid-state lasers), cam- 
eras (now holding two images recor ded in rapid su ccession) 
and computers (for post-processing) ( Adriani 120041) . The stan- 
dard system now usually comprises a single fast camera, a 
double pulsed Nd:Yag laser and a fast hard-wired PIV cor- 
relator. 

The basic method used to extract the velocity field from two 
successively acquired pictures can be decomposed in sev- 
eral steps. First, the raw image pairs are divided into several 
windows, called interrogations areas. Then, each pair of cor- 
responding windows ar e cross-correlated using F a st Fou rier 
Transform algorithms JScarano and RiethmuUen, Il999l) . A 
high cross-correlation value is observed where many par- 
ticles match up with their corresponding spatially shifted 
partners, and small cross-correlation peaks may be observed 
when individual particles match up with other particles. Only 
the highest correlation peak is interesting as its position in 
the correlation plane directly corresponds to the average par- 
ticle displacement within the interrogation area investigated. 
Finally, the velocity can be determinated, since the time sep- 
aration between the two pictures of one image pair is known. 
It appears that many parameters influence the velocity mea- 



surement and calculation: the number of particles image pairs 
in the interrogation area (lAdrianL Il99lh . the size of the par- 
ticles with respect to their spatial discretisation and the peak 
interpolation scheme. Recently, many efforts have been made 
to develop and improve more robust PIV algorithms: the 
spatial resolution and the measurement precision have been 
increased. The primary step was to reduce the size of the in- 
terrogation area from 64 x 64 p ixels or larg er ( Adriani 1 1 99 ll) 



to less than 10x10 pixels (.Nogueira et all 1200 11) and to use 



window offset to increase the spatial resolution (iWesterweel et al 



1997h . A secondary step was to reduce the uncertainty as- 
sociated with the interrogation procedure by using iterative 
methods (iScarano and RiethmuUen. Il999h . for example. As 
velocity measurements are necessary and now widely used 
to study turbulent flows, improvements a re con tinuously car- 
ried out to increase the spatial res olution ( Lavoie et al., 2007) 

or to reduce statistical errors (see Angele and Muhammad-KUngmann 
(l2005l) : ICholemaril ( l2007l) : lNogueira et alJhoOlh among oth- 
ers). 

However, if a higher spatial resolution can be needed to get 
a finer comprehension of turbulent flows, a higher time res- 
olution is also strongly necessary. For a long time, due to 
technical limitations, it was not possible to carry out high 
frequency PIV measurements (HFPIV). But recent improve- 
ments in the performances of computers, lasers and cameras 
make this development now possible. Recently, PIV systems 
operating at an acquisition rate up to 20 kHz have been u sed 
to study the wake in bluff bodies f Williams et al.L 12003^ or 
to obt ain velocity fields for both cold a nd hot flows (We rnet, 
20071). However the PIV system used by lWilliams et alJ (120031) 
consists in a copper vapour laser and a dump camera. This 
induces a heavy post-processing work and seems less adapted 
to the acquisition and the analysis of large series of velocity 
fields. 

In the present work, the flow field of a swirl-stabilized com- 
bustor, representative of Gas-Turbine combustor, is studied 
with a HFPIV system at an acquisition rate of 12 kHz. The 
structure of this paper is as follows. In section |2] the ex- 
perimental burner and the PIV system are described. Then, 
the structures of the mean reactive flow are detailed in sec- 
tion [5] so that the behavior of the burner can be studied and 
the efficiency and the resolution of the new system can be 
tested. Finally, instantaneous flow fields are presented and 
their analysis permits to gain deeper insight into the dynam- 
ical interactions between the flame, the vortices and the tur- 
bulent surrounding flow in the combustion chamber. 

2 Experimental Configuration 

2.1 Experimental Configuration 

Figure[T]shows a schematic diagram of the experimental in- 
jector. The setup is composed of a two-staged swirled in- 



jector, representative of a gas turbine injector, supplied with 
propane and air and a rectangular combustion chamber. 



the PIV system are reported on figure |2l 
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Fig. 1 Schematic view of the experimental injector. 

The primary stage consists in a central duct fed with pure 
fuel and a swirler supplied with air. The angle of the pri- 
mary stage swirler, which has 18 vanes, is maintained at 
42". The fuel/air mixture is realized in a mixing area down- 
stream the injections of propane and air. Fuel is delivered in 
the secondary stage by 15 holes. These 15 jets are located on 
a circular hollow part fed with propane. The airflow of the 
secondary stage is injected through a swirler with 20 vanes. 
The angle of the secondary stage swirler is 35°. Both air 
and combustible are mixed together in the secondary stage. 
Both swirlers turn the flow in the same direction, in order 
to obtain a strong co-rotating swirl motion. The primary and 
secondary fuel/air mixtures mix together while entering the 
combustion chamber. 

Both air swirlers are permanently supplied with air and the 
airflow rate in the secondary stage is 4 times higher than 
the first stage one. For the present study, for ease of illustra- 
tion, only the secondary stage fuel injection is used. Dried 
compressed air is available at a pressure of 0.7 MPa, while 
propane is stored at 0.5 MPa. Air and propane mass flows are 
monitored with electronic mass flow meters and controllers 
(Bronkhorst-Elflow). 

The combustion chamber has a square cross section of 100 x 
100 mm^ and a length of 500 mm. The side walls of the 
chamber are made of two silica windows in order to enable 
optical diagnostics in the flame, while the top and bottom 
walls are made of concrete. All tests are performed at ambi- 
ent pressure and ambient temperature. 



2.2 HFPIV measurements 

A High Frequency Particle Image Velocimetry (HFPIV) sys- 
tem is used to measure velocity flow fields in the combustion 
chamber. A schematic view of the experimental burner and 
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Fig. 2 Schematic view of the experimental combustor and the PIV sys- 
tem. 



For such experiments, the top and bottom walls of the cham- 
ber contain rectangular silica windows (25 mm wide and 
120 mm long) that allow the laser sheet to cross vertically 
the combustion chamber, downstream the injection plane. 
The light sheet is generated by a system consisting of two 
Nd: YAG lasers (Quantronix), a laser beam recombining de- 
vice and a set of cylindrical lenses (Melles Griot). Both lasers 
emit a pulse with a wavelength of 532 nm. They have a pulse 
energy and temporal width of 6 mJ and 160 ns respectively. 
Optics are used to combine both beams along the same tra- 
jectory. A set of lenses is used to transform the laser beam 
into a planar light sheet 60 mm wide and 0.3 mm thick. A 
fast speed camera (Photron Fastcam, 1024 x 1024 pixels at 
a rate of 2000 frames per second) equipped with a 105 mm 
F/1.8 Nikon Nikkor objective is placed perpendicularly to 
the burner axis to record the image pairs. 
The two lasers working at 12 kHz and the camera operating 
at 24 kHz are synchronised by a pulse delay generator {BNC 
555 pulses/delay Generator), as shown in the time diagram 
of figure [3] 
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Fig. 3 Time diagram for the time resolved PIV. 



At this rate of 24 kHz (for the camera), the size of the pic- 
tures is 128 X 512 pixels, a limited size which will be im- 
proved using new generation cameras. The airflow rate is 
seeded with Ti02 particles with a nominal diameter d = 
1 /xm. Ti02 has been used because its melting point is higher 
than the adiabatic flame temperature of the flame. These 
seeding particles will be present both in fresh and burnt gases 
Using estimations of the maximum velocity, the time delay 
between the two pulses is chosen equal to At = 6 /is. 
The raw image pairs are then exported and an off-line image 
processing is performed with an adaptative cross-correlation 
program ("Flow-Manager" by Dantec) using Fast Fourier 
Transform algorithm. The raw image pairs are divided into 
squared interrogation areas whose final dimensions are 8 x 8 
pixels with an overlap of 25 %. Finally, the spatial resolution 
of the flow field measured is 1 .07 x 1 .07 mm^. The PIV pro- 
gram is designed and optimized to deal with high velocity 
and density gradients, providing a peak finding error of less 



fer respectively to the horizontal and vertical velocity com- 
ponents. 



than 0.1 pixel in classical applications ( Westerweel et al 



19971). The parameters used for the HFPIV treatments are 
indicated in table [T| 



Table 1 PIV acquisition and processing parameters. All the images 
have the same size, however the region observed in the combustion 
chamber can be modified as described in section [3TT] D is the outer 
diameter of the injector 



Acquisition 


Field-of-view 


x = 0.47£);y=1.6fl 


Image size 


[128 512] pixels 


Camera pitch 


0.017 mm/px 


Frequency 


12 kHz 


Nb of raw images 


24,000 


Processing 


Sub-pixel scheme 


Gaussian 


Initial window size 


[32, 32] pixels 


Final window size 


[8, 8] pixels 


Overlap 


25% 


Nb of vectors 


[21, 85] 



2.3 Post-processing 

For each experiment, a continuous acquisition ofN^ 12000 
image pairs (24000 raw images) corresponding to a physical 
time of 1 second has been performed. The instantaneous ve- 
locity fields associated with the A^ image pairs are processed 
to determine the mean velocity components Umean and Vmean 
and the RMS velocity components Umis and Vnns (for "root 
mean square"). In the present study, the notation u and v re- 



i=\ 



^rms — \ \T 2^\^^ ^mean) 



(1) 



(2) 



!=1 



Each instantaneous velocity field is also processed in order 
to obtain the instantaneous vorticity field in the plane (x,y) 
using the following definition of the local vorticity, (O: 

du dv ^3^ 



(O 



dy dx 



where x refers to the horizontal direction and y to the verti- 
cal direction. As mentioned before, the spatial resolution of 
the velocity fields measured is rather low, making the cal- 
culation of du/dy and dv/dx difficult. Nevertheless, as it 
will be discussed in section HI this grid is fine enough to 
obtain interesting results concerning the behavior of the ex- 
perimental burner. 

While statistical quantities are necessary to get a good com- 
prehension of the burner behavior, time resolved results are 
also necessary to improve the understanding of the strongly 
unsteady phenomena occuring in the combustor As men- 
tioned in the introduction, these HFPIV measurements are 
mainly carried out to investigate the development of coher- 
ent structures in the combustion chamber. As, a clea r and 

universal definition of coherent strutures does not exist. ljeong and Hussain 
d 1995b consider the following to be the requirements for a 
vortex core: 

- a center of a vortex corresponds to a local pressure min- 
imum 

- a vortex core must have a net vorticity 

- the geometry of the identified vortex core should be Galilean 
invariant. 

Unfortunately, these three requirements do not permit to build 
a single identification scheme that can be used to analyse 
numerical or experimental data. It is usually difficult to ob- 
tain a pressure map experimentally, PIV measurements giv- 
ing only acces to velocity fields. A natural choice to extract 
coherent structures from PIV measurements would be the 
study of the instantaneous vorticity field, whose peaks coin- 
cide with coherent structures. However, it has been shown 

that vorticity cannot permit to distinguish bet ween swirling 

motion of a vortex from pure shearing m otions dJeong and Hussain , 
1995llKolaiil2007l:ISchram et al.ll2004 . This is why several 
methodologies based on the analysis of the velocity gradient 
tensor have been developed during the last t wo decades. A 
recent review of these methods is proposed in jKolan, 120071) . 
among others. 
One of the most widely used criterion to extract coherent 



structures from PIV meas urements, known as the X2 crite- 
rion, has been proposed by Jeong and HussainI d 1995b . They 
define a vortex in terms of eigenvalues of the symmetric ten- 
sor S^ + Q^ where S and Q are the symmetric and antisym- 
metric parts of the velocity gradient tensor Vu. As S^ + Q^ 
is symmetric, it has only real eigenvalues ordered as follows 
A 1 > A2 > A3. According to this criterion, a local minimum 
of the A2 field is associated with the local pressure mini- 
mum that is met in the vortex core. The study of the A2 field 
permits to discriminate between swirling and shearing mo- 
tions. However the calculation of the A2 field requires the 
velocity field and its variation in the three directions (x,y,z) 
which are not trivial to obtain experimentally. With 2D-P1V, 
only the two components of the velocity fields contained in 
the laser sheet plane and their variations in this plane are 
available. Because of this limitation and assuming that the 
flow is locally two-dimensional, which is a strong assump- 
tion here, the simplific ation of the A2 crite rion can be used 
JAnthoine et al., 2003:iSchram et al.ll2004 : 



(4) 



This criterion indicates that within the vortex, the derivatives 
of the velocity components must be of opposite signs. Thus, 
the diameter Dy of the vortex is defined by the zero crossing 
contour of A2. FinaUy, it has to be noted that the A2 field is 
a Galilean invariant, so that no additional transport should 
affect it. 

The instantaneous A 2 fields are first computed from the in- 
stantaneous velocity vector fields and then processed to ex- 
tract vortices. The detection method is based on the selectiv- 
ity property in space and scale of the wavelet transform of 
the in stantaneous A2 field (lAnthoine et al.Ll2003l : ISchram et al 
12004 ^. Assuming that the vortex signature in a A2 field is a 
Gaussian curve, this allows the determination of the vortex 
position and size. To detect such a shape, the Maars mother 
-wavelet has been selected to perform the wavelet analysis. 
The mother wavelet is given by: 



f(jt;,y) = (2-jt;2-/)exp 



-y 



(5) 



The wavelet transform of the A2 field is obtained by the con- 
volution product: 



W = (f,,,v|A2) = 



A2(A:,y)f/y yt/xc/y (6) 



where / is the scale of the mother wavelet, and a:' and y' 
are defined by: x' — x/l and y' = y/l. The detection algo- 
rithm first permits to calculate the wavelet coefficient, |W|, 
on each location (x,y) of the PIV plane. Then a vortex is 
identified if |W| is higher than a chosen threshold value and 
if it does not overlap another vortex. The detection algorithm 
provides the center of the vortices that are detected, assum- 
ing that each vortex center coincides with a local maximum 



of |W|. This makes it possible to estimate the velocity distri- 
bution around the vortices and their convection speed. The 
convection (or transport) velocity of a structure (m^; \\) is de- 
termined by averaging the velocity vectors contained in the 
structure, while the velocity distribution in the vortex core is 
obtained by substracting the convection velocity vectors to 
the instantenous field. These informations are necessary to 
follow vortex displacements in the combustion chamber and 
help in determining thek appearance frequency in the cham- 
ber. More details about the wavelet analysis can be found in 
refs JAnflioine et all 120031: Ischram et al.[l2004i) . 



3 Experimental results and interpretations 

3.1 Operating conditions 

The HFPIV diagnostic is tested for the operating conditions 
reported in table |2] 



Table 2 Operating Conditions. (Qc3H8)prim and (Qc3H8)sec are the 
fuel flow rates injected respectively in the primary and the secondary 
stages, while Qaj^ is the global air flow rate. <i> is the global equivalence 
ratio and P is the power delivered by the burner. 



Qair 


(Qc3W8)prim 


{Q_am)sec 


P 


4> 


Nm^.h-i 


Nm^.h-' 


Nm^.h-i 


kW 




105 





2.95 


74 


0.67 



To ease the analysis proposed in the present paper, only the 
second stage of the injector has been fed with propane. A 
complete analysis of the influence of staging on the burner 
stability will be proposed in a near future. 
The power associated with this regime is 74 kW. It may ap- 
pear low compared to industrial applications, but is already 
quite high for a laboratory scale burner A picture of the 
flame is shown on figure |4]for this particular regime. 




Fig. 4 Natural emission of the flame. P = 74 kW, g,„,- ^ 
Qcm» = 2.95 Nm\h-' and * = 0.67. 



105Nm\h- 



An experimental study of the acoustic behavior of the burner 
has shown that in the operating conditions chosen here, the 
burner exhibits strong combustion instabilities. This study 
has been carried out using a calibrated 1/4" microphone and 



a photomultiplier tube (PMT) equipped with a filter centered 
on the chemiluminescence band of CH* free radical (416 



nm) (IHiggins et alil200lh . The microphone (Bruel & Kjaer 
4136) placed in the bottom wall of the chamber measures the 
pressure fluctuations (p') on the semi-length of the chamber, 
while the PMT (Electron-Tubes, type 9124QB) collects the 
ligth emitted by the whole flame. The PSD of p' and the fluc- 
tuating part of CH* emissions [q') are reported on figure |5] 
Both PSD amplitudes present a strong peak, which corre- 
sponds to a frequency / == 316 Hz, meaning that thermo- 
acoustic instabilities are probably observed. 
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Fig. 5 PSD amplitude of microphone (dashed line) [dB] and PMT 
(Solid line) [logarithmic scale]. P = 74 kW, g„,v = 105 Nmlh"', 
Qam = 2.95 Nm^.h"' and (p = 0.67. 



This frequency can be associated with the quarter wave mode 
of the combustion chamber Indeed assuming that this cham- 
ber is homogeneously filled with burnt gases of average tem- 
perature T[,, the acoustic modes of the combustion chamber 
can be estimated as: / = (2n+ l)ci,/{4L), where Cj, is the 
sound velocity corresponding to 7], (cb = y/jrTi,), L is the 
combustion chamber length and n is the order of the reso- 
nant mode. It may be assumed that the average properties of 
the gaseous medium are equivalent to the ones of air. If we 
finally assume that lOOOK <Tb < 1500K, the quarter wave 
mode ( defined by n = 0) can be associated with frequencies 
between 300 and 375 Hz, a range which is consistent with 
the experimental results. It is expected from PIV measure- 
ments to gain insight into the dynamical coupling leading to 
combustion instabilities. 

Once the flow reaches a steady state and the flame is well 
stabilized in the combustion chamber, PIV images are ac- 
quired in the combustion chamber downstream the injection 
plane. Two complementary regions are defined for the mea- 
surements, as visible on figure |6] near the injection plane. 
The notations x* and y* are defined by x* ~ x/D and y* ~ 
y/D. The first region corresponds to an observation area 
which is 512 pixel high and 128 pixel wide and permits to 



visualize the flow field in a vertical window (0 < x* < 0.4 
and —0.7 <y* < 0.9). The second region is a 128 pixel high 
and 512 pixel wide window observation, that focuses on the 
phenomena that take place in the top part of the combustion 
chamber, between < .x:* < 1.8 and 0.4 < / < 0.87. While 
the first region makes it possible to analyze the distribution 
of vortices between the top and bottom parts of the cham- 
ber, the second one is devoted to the convection of structures 
along the chamber wall. 







Region 2 



1 Chamber 



I Axis 
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Fig. 6 Scheme of the two regions defined for the PIV measurements, 
with X* = x/D and y* = y/D. x* = corresponds to the begining of the 
combustion chamber, y* =0 coincides with the axis of the combustion 
chamber. 



Several validation routines have been used to perform data 
analysis: 

- Peak-height validation validates or rejects individual vec- 
tors based on the value of the peak height in the corre- 
lation plane where the vect or displacement is measur ed. 
The detectability criterion JKeane and Adrianl 119921) . k, 
that compares the highest validation peak with the sec- 
ond one is used to perform this validation method. 

- Velocity range validation rejects vectors which are out- 
side a certain expected range of velocities in the flow. 
This can be checked component by component or di- 
rectly on the velocity norm. 

- Moving-average validation validates or rejects vectors 
using a comparison between the vector tested and the av- 
erage of the vectors in a rectangular neighbourhood, us- 
ing the acceptance factor, a. This validation method is a 



partic ular case of iterative filtered validations (IHost-Madsen and McCluskey 
I994I) . 



As a result of a parametric study of the influence of va lida 



tion methods on the PIV results (iBarbosa et all 120081) . the 



values k—\.2 and 0: = 0.1 seems to be the most adapted to 
the present measurements. 

Nevertheless, results presented in the foHowing sections de- 
pend on the data post-processing and particularly on the val- 
idation methods used. The parametric study has shown that 
changing the validation parameters value introduces a max- 



imal uncertainty of 1 1 % in the maximal velocity, while it 
has no effect on the position of the maximal and minimal 
velocity regions dBarbosa et al.i , i2008i) . 



3.2 Mean flow field 

The mean axial and radial velocities, measured in a plane 
placed on the chamber axis, are reported on figure]?] The top 
of this vertical observation area is at the distance y* = 0.87, 
while the top (resp. the bottom) of the chamber is at a dis- 
tance of y* = 1 (resp. y* = — 1) and the edges of the injector 
are situated at y* = 0.5 (resp. — 0.5). 



Axial velocity |m/s] 



Radial Velocity [m/s] 




Fig. 7 Contours of the mean axial velocity U„,ean and mean radial ve- 
locity V,„j,„„. The velocity field is measured on the axis of the com- 
bustion chamber. The edges of the chamber and the injection area are 
drawn respectively in black and in gray. 



First, a large negative velocity area near the burner axis is 
observed on figure [T] It corresponds to the central recircula- 
tion region, which provides the major mechanism for flame 
stabilization in a swirl-stabilized burner. Second, the flow 
exhibits the features of an annular conical jet whose inter- 
nal boundary is visible on both sides of the cut at the dis- 
tances y* = 0.4 and —0.4. This observation indicates that 
the flame is clearly stabilized within the outer divergent of 
the injector The maxima of the mean axial and radial veloc- 
ities are observed within this annular jet. It can be noted that 
the maximum level of the radial mean velocity is almost one 
third of the axial component value. Two recirculation areas 
can also be observed in the top and the bottom corners of 
the chamber due to the sudden expansion at the injector exit 



plane. The annular jet induces two shear layers, which co- 
incide with the locations of the regions of maximum RMS 
velocity components, as visible on figure[8] As expected, the 
lower levels of turbulent intensities are obtained in the inner 
and outer recirculation areas. 

Axial RIvIS Velocity [m/s] Radial RMS Velocit,' [m/s] 
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Fig. 8 Contours of the axial RMS velocity J/,-,,,, and radial RMS ve- 
locity Vrms- The edges of the chamber and the injection area are respec- 
tively drawn in black and in gray. 



The distribution of the mean and RMS velocities are re- 
ported on figure |9] (axial component) and figure [TO] (radial 
component) for four successive axial cuts from x* =0.1 to 
X* = 0.4. 

The Umean profile, plotted in figure |9l is symmetric around 
y* = 0. At X* = 0. 1, it presents two peaks at a distance y* = 
±0.5, corresponding to the annular jet described before. These 
data will be used in section [33] to calculate the associated 
massflow rate and compare it with the value given by the 
massflow controller 

A large negative velocity area, corresponding to the inner 
backflow region, extends between y* — —0.35 andy* =0.35. 
Increasing x*, the location of the maximum level of Umean 
moves away from the axis of the injector, while the maxi- 
mum of Umean dccrcascs from an initial value of 44 m.s^' 
(x* = 0.1) to a final value of 37 m.s"' when x* = 0.4. Fig- 
ure|9]also permits to compare the distributions of {/,„«„„ and 
Urms- Profiles of Urms, plotted in figure |9] are symmetric 
around y* = 0. Urms presents maxima on each side of the 
chamber aty* ss ±0.43, i.e where the mixing region between 
the central recirculation zone and the positive streamwise 
velocity flow is located. Two secondary peaks, which are as- 
sociated with a far lower value and reached at y* « ±0.55, 
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Fig. 9 Profiles of the mean axial velocity Umean and the axial RMS 
velocity [/„„, for four successive vertical cuts, x* = 0.1, -t* = 0.2, x* = 
0.3 and X* =0.4. 



indicate the location of the external shear layer. When x* > 
0.2, the Urms profile presents only one peak in the mixing 
region between the recirculation area and the annular jet. In- 
deed, the outer recirculation has almost disappeared when 
X* > 0.2 (figureEJ. 

The same phenomena are visible on figure [TO] where the 
distributions of radial velocity Vmean and Vrms are reported. 
When X* < 0.2, the profile of Vrms presents two peaks dis- 
tributed on both sides of the annular jet as in the cse of Urms- 
However, whenx* > 0.2, Urms and Vrms do not follow similar 
trends anymore: the magnitude of Vrms increases monoton- 
icly with the radial distance (y* > 0), so that the peak of 
Vrms is observed near the top (resp. the bottom) of the com- 
bustion chamber Due to the coarse spatial resolution of the 
measurement technique, this point needs to be further inves- 
tigated. 
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Fig. 10 Profiles of the mean radial velocity V„,ja„ and the RMS of ra- 
dial velocity Vn„s for four successive vertical cuts, x* = 0.1, x* = 0.2, 
X* =0.3 and X* =0.4. 



In order to better visualize the phenomena occurring around 
the jet, the camera focuses on the top of the chamber (0 < 
^* < 1.8 and 0.4 <y* < 0.87). Figure[TT]displays the mean 
velocity field obtained by processing A^ = 12000 instanta- 
neous fields. The central back-flow region, revealed by the 
streamline plot and the negative axial velocity area, extends 
between x* = 0.4 and 1.8, while the second recirculation 
area, in the top of the chamber (y* > 0.6) develops between 
X* = and 0.3. The annular jet characterized by positive ve- 
locity vectors separes these two recirculation areas. 
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Fig. 11 Streamlines of the mean velocity field are plotted on the top 
part of the figure. Only the top part of the combustion chamber is repre- 
sented. The mean axial and radial velocity are displayed in the middle 
and the bottom parts of the figure. The velocity field is measured on 
the axis of the combustion chamber. 



The distributions of the mean axial and radial velocities ob- 
tained with a vertical observation area ("Region 1") and a 
horizontal one ("Region 2") are plotted on figure [T2]for the 
vertical cut x* =0.1. The two measurements show a very 
good agreement in terms of Umean and Vmem- The observa- 
tion area has no effect on the position of the velocity peaks, 
while a difference in amplitude of less than 5 % is observed 
in the region of high velocities, which is quite acceptable. 
Using this recording, the evolution of the maximum of Unean 
in each vertical cut is plotted in figure[T3]as a function of the 
reduced axial distance to the injection plane. On the same 
graph, the evolution of the radial position of this maximum 
velocity is reported using the left ordinate axis. 
The Umean peak level increases initially up to a distance of 
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Fig. 12 Profiles of V„,f,„„ (rigth) and U„„„„ (left) for;i:* =0.1. Profiles 
obtained with the vertical observation area ("region 1") are plotted in 
gray and the ones obtained with the horizontal observation area ("re- 
gion 2") are plotted in black. 




X* — 0.2. The profile of Umis downstream the injection plane 
present two peaks at x* =0.1 and 0.25, corresponding to the 
location of the external and internal shear layers due to the 
annular conical jet. Close to the wall, Urms presents a peak, 
which seems to have no physical signification and might be 
due to the boundary effects. The profile of [/,„„ indicates 
that the maximum of fluctuations takes place downstream 
the injection plane: x* < 0.6. Farther downstream (x* > 0.6), 
the mean axial velocity becomes negative and Urms reaches 
its lowest level. This indicates that the internal recirculation 
area extends almost beyond 2 diameters from the burner. 



3.3 Calculation of the flow rate injected 

In order to check the coherence of the experimental results, 
the fluid flow rate in the combustion chamber is calculated 
using the axial velocity measured experimentally using HF- 
PIV. The fluid flow rate is estimated using the radial distri- 
bution of Umean, taken at a distance of x* = 0.06, so that 
the result is less disturbed by edge effects. Assuming that 
the flow is symmetric, the flow rate is calculated using the 
following definition: 



Fig. 13 Evolution of the maximum value of £/,„„,„ with x* . The evolu- 
tion of the normalized radial position corresponding to this maximum 2 = / / V.nc/5 



X* — 0.2 where it reaches 43 m.s^' for a transverse distance 
y* ~ 0.5. Further downstream the U,„ecm level decreases from 
43m.s^' to 25m.s^', while the radial position of the peak 
velocity increases fromy* = 0.5 to 0.8. It can be considered 
as indicative of the position of the flame front in the cham- 
ber, assuming that the positive streamwise velocity separates 
the flame from the fresh gases. 

The distributions of U,„ean and Urms are finally plotted on 
figure[T4]in the shear layer vicinity (y* = 0.5). 



(7) 



where v is the velocity vector Imposed and measured results 
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Fig. 14 Evolution of Umean and Umu withx* fory* = 0.5. 



As already mentioned, Umean reaches its maximum level at 



are respectively Qtk = 113 m .h^ and Qexp — 115 m.h 
meaning that the relative error between imposed and mea- 
sured flow rates is about 2% for the case studied. This dif- 
ference can be due to the precision of the mass flow me- 
ters, the precision of the PIV measurements, the assump- 
tion of axisymmetry of the flow or the position of the profile 
used. However the relative error is weak, and the experi- 
mental results seem consistent. Finally, this good agreement 
can be used to say that PIV measurements give a good de- 
scription of the phenomena taking place in the combustion 
chamber on a quantitative point of view. The mean flow field 
structures have been characterized using PIV measurements. 
However, the principal advantage of this diagnostic is the 
high acquisition rate of the image pairs, which are then pro- 
cessed to calculate instantaneous velocity fields. An analy- 
sis of the instantaneous flow field evolution in the chamber 
is reported the following section. 



4 Analysis of instantaneous flow fields 

A temporal analysis of the velocity field is carried out in the 
regime defined in tablelH A serie of six successive instanta- 
neous velocity vectors and their associated vorticity field are 
displayed in figure [15] The views are taken on the vertical 
symmetry axis of the combustion chamber (y* = 0). A time 
delay of At = 83 /is separates each recording. This serie has 
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Fig. 15 Instantaneous velocity vectors and associated vorticity fields for six successive recordings separated by a time delay of 83 /is. These fields 
are obtained on the axis of the combustion chamber. On each picture, the location of the injector outlet is represented in black. 



been chosen because it is representative of the phenomena 
observed. 

The instantaneous velocity fields exhibit large deviations from 
the mean velocity field described on figures [T] and [8] Still, 
a large inner recirculation area and an annular conical jet 
are clearly visible on each picture. The inner recirculation 
area, characterized by a negative velocity and a low vortic- 
ity level, extends on the center of each picture. The location 
of the boundary of the annular jet changes dramatically from 
one instantaneous image to the others. Figure [15] shows that 
maximum levels of vorticity, o), coincide with the jet bound- 
aries, where the flow is characterized by a great shear rate. 
Finally, a large scale vortex is visible at y* w 0.4 on the first 
frame of figure[T5] The five other frames show the deforma- 
tion and convection of this structure, and a modification of 
its vorticity value. 

As seen before a straightforward criterion to identify a vor- 
tex is to detect the presence of a vorticity peak. However, 
while the center of a vortex is actually associated with a 
vorticity maximum, a vorticity peak does not necessary cor- 
respond to a vortex center. Indeed, on figure [Tsl high val- 
ues are present on the external shear layer surrounding the 
annular jet, where no vortex can be detected. An analysis 
based on the A 2 method presented by Jeong and Hussain 
Jeong and HussainI ( Il995h and describe in paragraph 12.31 is 
then proposed to identify vortices. 



4. 1 Vortex detection 

Figure [T6l shows the A 2 fields, analyzed with the wavelet 
transform, for the six velocity fields of figure [15] The A 2 
criterion is not sensitive to shear regions, on the contrary to 
vorticity. Several peaks of A 2 are detected and the structure 
cited in the previous paragraph (figure [Tst is visible in the 
top part of the combustion chamber (y* w 0.4). 
Using the detection algorithm described in paragraph 12.31 
the axial displacement and the convection velocity of this 
structure are determined and their evolution as a function of 
time is reported on figure [TT] 

First, for / < 160/Xs, the axial convection speed Uc is large 
and positive, explaining the displacement of the structure 
towards the right end of the observation area. When, t > 
160 /is, the value of Uc becomes slightly negative explain- 
ing that the location of the structure does not change, as if 
the structure was trapped at the limit of a region where the 
flow recirculates. The value of the radial transport velocity, 
Vc, is almost constant, close to 12 m.s^'. 
In the 12000 instantaneous velocity fields of the test, the 
phenomena described on figures [15] [16] and [TT] have been 
periodically observed on both half parts of the chamber. The 
values obtained on figure[T7]give a good estimate of the con- 
vection speed of structures created within the combustion 
chamber. Lastly, as A2 is a Galilean invariant, the detection 
algorithm permits to find small vortices that are difficult to 
visuaUse due to the high levels of velocity encoutered in the 
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Fig. 16 Velocity vectors and associated wavelet transform of the 12 field for the six successive recordings shown on figure [Tsl For At — Qfls, the 
velocity field around each structure detected is shown. On each frame, the location of the injector outlet is represented in black. Left: Enlargement 
of two zones where vortices have been detected. 
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Fig. 17 Distribution of the axial position, the axial and radial transport 
velocities as a function of time, for the six fields of figure [16] The 
radial position of the followed structure is y* « 0.4. The time t = fis 
coincides with the first frame of figure [T6l 



chamber as shown on the details of the first frame of fig- 
ure [16] In spite of the coarse grid used to calculate the ve- 
locity vector field, these structures are detected. 



4.2 Coherent structure frequency 

The coherent structure frequency is another important pa- 
rameter for the burner characterization. As the detection al- 



gorithm gives the coordinates of each detected vortex, it is 
possible to count the number of vortices that crosses a given 
inteiTogation area of the chamber during a given time ( 1 sec- 
ond for the present study). For a given observation area, if 
the value of |A2| is higher than a given threshold a vortex 
is detected. The number of the frame and the value of |A2| 
is noted so that it is possible to plot the peak of |A2| versus 
time as visible on figure[T8] In order to illustrate phenomena 
taking place in the chamber, a squared observation area of 
3 mm wide is centred on the point of coordinates {x* = 0.24, 
/ = 0.42). 

On figure [18] all peaks of \k2\ have been normalized using 
the maximum value encountered during 1 second of record- 
ing. The analysis of this interrogation area indicates a vortex 
frequency of about 310 Hz. The same analysis has been per- 
formed using several observation area located on the top or 
the bottom part of the annular jet and the frequency of vor- 
tices has always been found between 300 and 320 Hz, which 
corresponds to an uncertainty of about 3 %. To complete this 
analysis, the power spectral density (PSD) of the instanta- 
neous axial and radial velocities are calculated in each point 
of the velocity grid. For each experiment, the acquisition 
frequency is / = 12000 Hz and the number of samples is 
A^ = 12000, so that the spectral resolution is Af = 1 Hz. 
To illustrate the results of the burner aerodynamic spectral 
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Fig. 19 PSD amplitude of the axial and radial velocities at point A of coordinates (x* = 0.2, y* = 0.4) and point B of coordinates (x* = 0.2, 
/ = -0.4). 
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Fig. 18 Vortices detected within a squared observation window 3 mm 
wide centred on the point of coordinates (x* = 0.24, y* = 0.42). 



analysis, the velocity PSD amplitudes calculated at the point 
A (jc* = 0.2; y* = 0.4) and the point B (x* = 0.2; y* = -0.4) 
are plotted versus frequency in figure[T9l The points A and B 
are located at the boundary between the central re-circulation 
area and the annular jet. As seen before vortices are mainly 
created in this region of great shear rate. These two points 
are chosen symmetric around the chamber axis. 
The PSD amplitude of u and v calculated at point A are re- 
ported on the left part of figure [19] and the same analysis has 



been performed in the bottom part of the combustion cham- 
ber, for point B. The PSD amplitudes of the axial and radial 
velocity calculated are reported on the right part of figure[T9l 
The PSD amplitude of u indicates that the burner behaves 
like a low-pass filter: the highest levels of the PSD amplitude 
are observed in the range of frequencies < / < 1000 Hz 
and then the PSD amplitude decreases to reach its lowest 
level when / > 2000 Hz. The u spectra for A and B are 
quite similar and they both present a peak at a frequency 
/„«310Hz. 

The main peak value of the v spectrum is also observed at a 
frequency /,, ss 310 Hz for point B, even if this peak is less 
well marked. This last point is even more pronounced for 
the V spectrum of point A, where the signal is broadband in 
the range [0; 1000] Hz 

It can be noted that /„ and /y are very close to the vortex fre- 
quency determined before and to the acoustic quater wave 
mode of the combustion chamber (paragraph [TT)- One can 
imagine a straightforward scenario, in which the flame heat 
release rate is pulsated by strong aerodynamic structures, 
whose frequency is the first eigenmode of the combustion 
chamber. 
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5 Concluding Remarks 



References 



The aim of the present study was to detail and evaluate the 
efficiency of a system of HFPIV operating at 12 kHz. This 
diagnostic is used to characterise the behavior of an experi- 
mental lean premixed swirl-stabilized burner representative 
of a gas turbine combustor, which may exhibit strong com- 
bustion instabilities under certain operating conditions. It is 
expected from HFPIV to provide well temporally resolved 
velocity fields that permit to better understand these com- 
plex turbulent reactive flows. 

Using HFPIV measurements, the mean structure of the flow 
field has been first determined. This mean velocity field ex- 
hibits the features of a swirling flow, with an annular con- 
ical jet where high values of the axial velocity component 
are observed, and a large central recirculation zone that pro- 
vides the major mechanism for flame stabilization. RMS ve- 
locities are also estimated, showing large levels of fluctua- 
tions inside the combution chamber, associated with a strong 
dynamical behavior of the flame especially in the shear re- 
gions. These mean quantities are also used to estimate the 
efficiency of the HFPIV system. The mass flow rate enter- 
ing the chamber is calculated and compared with the value 
given by the mass flow controllers and a very good precision 
is obtained (error around 5 %). 

Instantaneous velocity fields are analyzed in order to better 
understand the interactions between vortices and the turbu- 
lent flame in the combustion chamber. Large coherent struc- 
tures are periodically detected within the combustion cham- 
ber. Their convection speed and their frequency are deter- 
mined and it is shown that, in resonant situations, this fre- 
quency can be related to an acoustic resonant mode of the 
combustion chamber. 

To conclude, to achieve temporally resolved measurements, 
it has been necessary to find a compromise between spatial 
and temporal resolutions, due to the technical limitations of 
high acquisition rate cameras at the present time. The study 
shows that the quite low spatial resolution is high enough 
to characterize the flow and to detect coherent structures. 
Finally, the present work indicates that the high temporal 
resolution chosen (12 kHz) is necessary to study turbulent 
GT flames where strongly unsteady phenomena occur. The 
actual technical developments of cameras lead to suppose 
that more interesting compromises between image size and 
acquisition rate can be found in the next years, so that it 
will become possible to obtain temporally resolved velocity 
fields measurements with HFPIV, with a better resolution. 
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